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ABSTRACT 



Recent INTEGRAL observations and archival BeppoSAX 
data have been analyzed to deeply investigate the hard 
X-ray behavior of the neutron star, atoll type, low mass 
X-ray Binary 4U 1636-53 . Our investigation in three dif- 
ferent periods outline three corrisponding different sates. 
Infact, the source was detected in both the canonical high 
and low state and moreover in one occasion INTEGRAL 
spectrum shows, for first time in this source, a hard tail 
dominating the emission above 30 keV. This spectrum is 
fitted as the sum of a Comptonized component similar to 
soft state and a power-law component (r = 2.76), indi- 
cating the presence of a non thermal electron distribution 
of velocities. 
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1. INTRODUCTION 



4U 1636-53 is a neutron star low mass X-ray binary 
(LMXB) classified as a atoll source IU7I1 . with an or- 
bital period of 3.8 hr derived from the optical variabil- 
ity of its companion V801 Arae l25ll and at distance of 
3.7-6.5 kpc 1151 1271 y]. While the X-ray burst proper- 
ties and timing signatures have been analyzed extensively 
(see lUsTl . 0] and references therein) the spectral charac- 
teristics have been studied only at low energy with Ein- 
stein it, EXOSAT HI, Temna 01 and ASCA 0. In 
general, the spectrum was acceptably fitted by a Comp- 
tonization model plus a black body component. We re- 
port here a broad band spectral analysis performed on 
data from BeppoSAX and INTEGRAL satellites, which 
allowed us to better constrain the spectral parameters 
and to detect the presence of a high energy tail domi- 
nating the spectrum above ~ 30 keV. A similar feature 
has been observed in other LMXBs namely GX 17+2 
lO, GX 349+2 0, Sco X-l U, 4U 1608-522 El, 
XB 1254-690 |20], Cir X-l [21] and 4U 0614+091 [26]. 



2. OBSERVATIONS AND DATA ANALYSIS 



TableHsummarizes the log of INTEGRAL and BeppoSAX 
observations of 4U 1636-53 . BeppoSAX observed the 
source on three occasions: February and March 1998 
and February 2000. LECS, MECS and PDS event files 
and spectra, available from the ASI Scientific Data Cen- 
ter, were generated by means of the Supervised Standard 
Science Analysis J3. Both LECS and MECS spectra 
were accumulated in circular regions of 8' radius. The 
PDS spectra were obtained with the background rejec- 
tion method based on fixed rise time thresholds. Publicly 
available matrices were used for all the instruments. Fits 
are performed in the following energy band: 0.5-3.5 keV 
for LECS, 1.5-10.0 keV for MECS and 15-70 keV for 
PDS. 

The analyzed INTEGRAL | 33] data set consists of all ob- 
servations in which 4U 1636-53 was within the high- 
energy detectors field of view. Observations are orga- 
nized into uninterrupted 2000 s long science pointing, 
windows (sew): light curves and spectra are extracted 
for each individual sew. Wideband spectra (from 5 to 
150 keV) of the source are obtained using data from the 
two high-energy instruments JEM-X 12311 and IBIS 1 29]. 
Data were processed using the Off-line Scientific Anal- 
ysis (OS A version 5.1) software released by the INTE- 
GRAL Scientific Data Centre. Data from the Fully Coded 
field of view only for both instrument have been used. 



3. SPECTRAL ANALYSIS 



Several physical models have been used to fit the whole 
data set, keeping the number of free parameters as low 
as possible. Each time a new component was added to 
the model, a F-test was performed. We assumed that a 
F probability larger than 95% implies a significative im- 
provement of the fit. The uncertainties are at 90% confi- 
dence level for one parameter of interest (A% 2 = 2.71). 
When spectra are from more than one detector, we al- 
low the relative normalization to be free with respect to 
the MECS and IBIS data, for BeppoSAX and INTEGRAL 
respectively. XSPEC v. 11.3.1. has been used. 



Table 1. BeppoSAX and INTEGRAL Observations 
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Spectral behavior has been studied separately in three 
epochs consisting of the following data: 

• I s * epoch: all three BeppoSAX observations avail- 
able from February 1998 to February 2000. During 
these periods the source was always in a soft/high 
state. 

• 2 nd epoch: JEM-X and IBIS data available from 
52644 MJD to 53644 MJD with count rate < 
5 counts s _1 in the 20^-0 keV energy band. For 
the chosen period the source was in a hard/low state. 

• 3 rd epoch: JEM-X and IBIS data available from 
52644 MJD to 53644 MJD with count rate > 
5 counts s _1 in the 20^0 keV energy band. This 
epoch does not corrspond to either the soft or the 
hard state and and here we call it peculiar state as 
will be explained in detail later on. 

After veryfing there were no changes in the three obser- 
vations and in order to achieve the highest signal to noise 
ratio, we build the BeppoSAX average spectrum arrang- 
ing all three observations. We then can take advantage 
of the high quality of the average spectrum up to about 
70 keV. The most simple model which provides a good 
fit to the average BeppoSAX spectrum in the energy band 
0.5-70 keV consist of a thermal Comptonized component 
modeled in XSPEC by comptt (a spherical geometry was 
assumed) plus a soft component which we modeled with 
two temperature blackbody. Being the column density 
Nh towards the source measured by the LECS and MECS 
instrument always close to the galactic column density 
(Ah galactic = 3.58 x 10 21 cm -2 , estimated from the 
21 cm measurement 1 10]), we fix this parameter at this 
value. Spectral fit results are given in Table the aver- 
age spectrum is shown in Fig. 1. 

The source was in the soft/high state with an un-absorbed 
luminosity of £0.1-200 keV — 2.0 x 10 37 erg s -1 , as- 
suming a distance of 5.9 kpc (Cornelisse et al. 2003). 
As in the case of the Einstein observation [6], thermal 
Comptonization of the optically thick plasma corona with 
a quite low electron temperature is dominating. 
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Figure 1. The BeppoSAX average spectrum in the soft 
state ( epochs 1 ), shown together with the total model and 
its components: the total, two blackbody and the Comp- 
tonization components are shown in green, red, magenta 
and blue, respectively. 

The most simple model which provides a good fit to 
the INTEGRAL hard state consist of a thermal Comp- 
tonized component modeled in XSPEC by comptt l28Tl 
(a spherical geometry was assumed) plus a soft compo- 
nent which we modeled by a single temperature black- 
body. Because of the very good low energy BeppoSAX 
coverage, we used the average value of input soft pho- 
ton temperature and column density (To = 1.3keV and 
Nh = Nh Galactic) measured by the BeppoSAX ob- 
servations. During this period, the source has a luminos- 
ity of Lo. 1-200 keV — 1-4 x 10 37 erg s _1 , lower than the 
one in the soft state, in agreement with the usual ranking 
of the luminosity in atoll sources (e.g., lTl7ll : ll3l1l : flfilo . 
The electron temperature is now substantially higher than 
in the soft state, T c ~ 23 keV, and the Comptonization 
component extends well above ~ 100 keV. The same 
model was been applied to the 2 nd INTEGRAL data set 
but resulted in a poor fit with a x 2 /d.o.f — 97/57 with 
clear residuals above 60 keV. Adding a power law im- 
proves the fit significantly (x 2 /d.o.f becomes 64/55). 
The disk component becomes negligible and it is not nec- 
essary to best fit the data. Spectral fit results are given 
in Table [2] and spectra are shown in Fig. |2] Since the 



Table 2. Results of the fit of 4U 1636-53 spectra in the energy band 0.5 — 70 keV and 5 — 150 keV for BeppoSAX and 
INTEGRAL observations respectively, for three different spectral states. 
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spectral parameters in the peculiar state are similar to 
those in the soft/high BeppoSAX state, a possible inter- 
pretation of this state is that 4U 1636-53 was in a simi- 
lar soft/high state as during the BeppoSAX observations, 
but with a new overlapped component at high energies 
simultaneously present. The un-absorbed luminosity is 
£o.i-20okcV ~ 7.4 x 10 37 erg s -1 . 



4. CONCLUSIONS 



The black-body component in the soft states could orig- 
inate at both the neutron star surface and the surface of 
an optically-thick accretion disk. In our observations the 
two black body components seems to originate from two 
different parts of the disk, corrisponding to two differ- 
ent temperature. The Comptonization component may 
arise from a corona above the disk and/or between the 
disk and the neutron star surface. In the hard states, ac- 
cretion probably assumes the form of a truncated outer 
accretion disk as previously reported by |3] for LMXB 
4U 1705-44. The spectral transitions are generally, but 
not necessarily, coupled with changes in luminosity, in- 
dicating they are driven by variability of the accretion 
rate or change of the geometry of the system as for Black 
Hole hard/soft transition at constant luminosity |4|. For 
4U 1636-53 the accretion rate is lower in the hard state 
than in the soft state and becames very high in the pecu- 
liar state. We note this value can be influenced by the 
steep power law that dominates the energy spectrum at 
low energies. The emission is best described by Comp- 
tonization from a complex electron distribution due to 
a low temperature (~ 5keV) thermal electron distribu- 
tion together with non-thermal power-law electrons. This 
two-component electron distribution could be explained 
by the following hypothesis: 

• Non-thermal electron acceleration regions powered 
by magnetic reconnections above a disc. Low- 
energy electrons cool preferentially by Coulomb 
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Figure 2. The spectra of epochs 2 ( soft state, left) and 
3 (peculiar state, right) observed by INTEGRAL , shown 
together with the total model and its components. Left: 
the total, the blackbody and the Comptonization compo- 
nents are shown in green, magenta and blue, respectively. 
Right: the total, the Comptonization and the power law 
component are shown in green, blue and orange respec- 
tively. 



collisions leading to a thermal distribution while the 
high-energy electrons cool by Compton scattering, 
preserving a non-thermal distribution [3]. 

• The magnetic reconnection above the disc can pro- 
duce a non-thermal electron distribution, while over- 
heating of the inner disc produces the thermal 
Comptonization 1 19]. 

• The power-law component could be produced by 
Componization by synchrotron emission in the jet 
0, loi In our source, this hypothesis is strength- 
ened by radio detection: Sydney University Molon- 
glo Sky Survey catalog gives a flux of 7.5 mJy at 
843 GHz |H. 

Finally, Laurent & Titarchuk 1 22] suggest the detection 
of a power-law component at high energy to be a signa- 
ture of presence black hole in an X-ray binary system. 
Our data are supporting result on GX 17+2 Hill , clearly 
showing this criterion proposed to distinguish black hole 
versus neutron star binaries is inadequate. 
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